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Magnetic knots as the origin of spikes in the gravitational wave backgrounds 
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The dynamical symmetries of hot and electrically neutral plasmas in a highly conducting medium 
suggest that, after the epoch of the electron-positron annihilation, magnetohydrodynamical config¬ 
urations carrying a net magnetic helicity can be present. The simultaneous conservation of the 
magnetic flux and helicity implies that the (divergence free) field lines will possess inhomogeneous 
knot structures acting as source “seeds” in the evolution equations of the scalar, vector and tensor 
fluctuations of the background geometry. We give explicit examples of magnetic knot configura¬ 
tions with finite energy and we compute the induced metric fluctuations. Since magnetic knots 
are (conformally) coupled to gravity via the vertex dictated by the equivalence principle, they can 
imprint spikes in the gravitational wave spectrum for frequencies compatible with the typical scale 
of the knot corresponding, in our examples, to a present frequency range of Hertz. At 

lower frequencies the spectrum is power-suppressed and well below the COBE limit. For smaller 
length scales (i.e. for larger frequencies) the spectrum is exponentially suppressed and then irrele¬ 
vant for the pulsar bounds. Depending upon the number of knots of the configuration, the typical 
amplitude of the gravitational wave logarithmic energy spectrum (in critical units) can be even four 
orders of magnitude larger than the usual flat (inflationary) energy spectrum generated thanks to 
the parametric amplification of the vacuum fluctuations. 
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I. INTRODUCTION 

The idea that large scale magnetic fields could exist in the early Universe is both appealing and plausible Q . It was 
noticed long ago that if magnetic fields existed in the early Universe, they could have affected (or possibly explained) 
various stages of the formation of cosmic structures Long range stochastic fields can modify the rate of the 

Universe expansion and the reaction rate at the nucleosynthesis epoch 1^. Stochastic fields coherent over the horizon 
size at the decoupling time can depolarise the Cosmic Microwave Background Radiation (CMBR) [^, and can also 
change the sound velocity of the plasma at the corresponding era [Q . 

If the temperature of the radiation-dominated Universe is higher than the temperature of the electroweak phase 
transition (i.e. T > 100 GeV), then the large scale components of the standard model gauge fields should be identified 
with the hypermagnetic fields which, by contributing to the (Abelian) anomaly, can generate matter-antimatter 
fluctuations ^ (possibly relevant in the context of the inhomogeneous nucleosynthesis scenario |^,|l^). Long range 
hypermagnetic fields can also influence the dynamics of the phase transition itself j^. Interesting consequences of 
the presence of electroweak condensates for the gravitational wave spectra were also discussed in . 

Stochastically distributed magnetic fields (either amplified from their vacuum fluctuations jl^ through the breaking 
of the conformal invariance of Maxwell equations or just assumed [T^) can imprint inhomogeneities in the microwave 
sky. Long ago it was argued that any uniform (spatially homogeneous) magnetic fields can make the energy momen¬ 
tum tensor of a radiation dominated universe slightly anisotropic producing, ultimately, an observable effect in the 
microwave sky (^. Recently bounds were derived (using the COBE 4-yrs maps) on the present strength of a uniform 
(and homogeneous) magnetic field ]T^ . 

In this paper we will argue that, on the basis of the physical laws governing the evolution of the mean magnetic 
fields in a hot plasma, there are no reasons why the topology of the magnetic flux lines should be trivial, namely 
there are no reasons why the magnetic flux lines should not intersect each others. If the topology of the magnetic flux 
lines is non trivial, then, the helicity of a specific magnetic configuration will also enter in the equations describing 
evolution of the scalar, vector and tensor fluctuations of the metric in a radiation dominated Universe. Provided the 
typical scale of the knot (red-shifted at the recombination time) is much smaller that the magnetic Jeans scale at the 
corresponding epoch, very little impact can be expected for structure formation ]l7| , p^ . 

In short our logic is the following. One of the main features of the evolution equations of the mean magnetic 
fields (for scales larger than the Debye radius) in highly conducting (and globally neutral) plasmas is that they are 
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divergence-less (i.e. the field lines have no end points) and, therefor^the structures arising in magnetohydrodynamics 
(MHD) can be examined in terms of the topology of closed curves Now, it is well known that in (ideal) MHD 

a unique and intriguing property holds: the magnetic flux and magnetic helicity are simultaneously conserved all 
along the evolution of the plasma for scales larger than the magnetic diffusivity scale. Since the Universe is a good 
conductor all along its history, we have to conclude that, on one hand, thanks to the magnetic flux conservation, 
the magnetic flux lines always evolve glued together with the plasma element and, on the other hand, thanks to the 
magnetic helicity conservation, the sum of the number of twists and the number of knots in the flux lines is also left 
invariant as time goes by. At this point the choice is twofold. Either we assume (as it is frequently done in the 
case of stochastic magnetic backgrounds @11) that the (mean) magnetic helicity is strictly zero, or we accept that 
the magnetic helicity can be present and, in this last case, it will be conserved to a very good approximation. 

As we said, a non-vanishing magnetic helicity implies that the corresponding magnetic field configurations must 
be inhomogeneous (otherwise the integral defining the magnetic helicity cannot be gauge-invariant). Thus, in the 
present paper, we are forced to deal with inhomogeneous fields. At the same time one could say that inhomogeneous 
configurations carrying zero (net) magnetic helicity can also have interesting effects on the gravitational wave produc¬ 
tion. This is certainly true and one of the purposes of the present investigation is to understand, with some specific 
examples, what difference does it make to drop the (oversimplifying) assumption that the magnetic flux lines of the 
cosmological magnetic fields have all trivial topological structure. 

Inhomogeneous magnetic knots are expected to affect various aspects of the life of our Universe. They can represent a 
quite interesting source “seed” in the evolution of metric density perturbations during radiation (and possibly matter) 
dominated epochs. Since magnetic knot configurations with finite energy and helicity are necessarily inhomogeneous, 
no homogeneous and isotropic (background) contribution is expected from the knot energy-momentum tensor. At the 
same time, the presence of magnetic knots and twists in the magnetic flux lines does modify the evolution equations of 
the metric fluctuations. The spatial components of the knot energy-momentum tensor form a three-dimensional rank 
two (cartesian) tensor which can be decomposed into irreducible scalar (invariant under three dimensional rotations), 
vector and tensor parts providing source terms for the evolution equations of the corresponding scalar vector and tensor 
parts of the metric inhomogeneities. The tensor modes of the metric are automatically invariant under infinitesimal 
coordinate (gauge) transformations pl[| . The components of the knot energy-momentum tensor are also invariant 
under infinitesimal coordinate transformations because of the absence of any (homogeneous and isotropic) magnetic 
background. 

We want to notice immediately that, one of the main ideas of our investigation is the connection among knot seeds 
and metric fluctuations. Up to now all the investigations on the existence of large scale magnetic fields were mainly 
devoted to the study of suitable mechanisms able to explain the large scale “magnetic” structure of the Universe, or 
more specifically, the existence of galactic (and possibly inter-galactic) magnetic fields @] . In our study we point out 
that, if the conservation properties arising in MHD are properly taken into account, then, interesting effects can be 
foreseen also for the metric fluctuations. 

In the first place we would like to motivate magnetic knots configurations on the basis of magnetohydrodynamics 
(MHD). The idea of magnetic knots might seem, at hrst sight, a crazy one. In order to make this idea quantitatively 
plausible we will introduce some explicit MHD solutions describing magnetic knots. It is very important that these 
field configurations are regular over the whole space, they have finite energy and helicity. They represent then a good 
framework where the effect of magnetic knots on the metric inhomogeneities can be discussed. The fact that our 
configurations have finite energy (as also required by considerations related to the C/em(l) gauge invariance of the 
magnetic helicity) implies that the amplitude of the magnetic fields tends to zero for scales larger than the scale of the 
configuration Ls- For the discussion of the magnetic helicity evolution in a radiation dominated stage of expansion 
an important tool will be the conformal invariance of the (ideal) MHD equations. Moreover, since we want our MHD 
description to be fully valid we will also consider times larger than the epoch of electron-positron annihilation. 

The second purpose of this investigation are the possible phenomenological consequences of the existence of magnetic 
knots prior to the decoupling time. We will show that spikes in the gravitational waves energy spectrum p^ ] can 
be foreseen and we will compare our results with the typical amplitudes of gravitational waves produced in the 
context of inflationary models. The amplitude of the “magnetic” spikes depends directly upon the helicity of our 
configuration and turns out to be substantially larger than the signal of a stochastic gravitational wave background with 
flat (Harrison-Zeldovich) logarithmic energy spectrum which represents a generic prediction of ordinary inflationary 
models. The “magnetic spikes” in the gravitational wave spectrum are not constrained by pulsars’s pulses since 
they are expected to occur in a (present) frequency range ujdec < oj < oj^+e- (where uJdec ~ 10“^® Hz and We+e- ~ 
10“^° Hz. In principle we could also expect interesting effects from the Lorentz force term associated with the knot 
configuration and appearing in the Euler (Navier-Stokes) equation for the density field 01] which goes, in MHD, 
as J X ~ (V X H) X H. The calculation shows that, since the knot scale is just slightly smaller than the magnetic 
Jeans scale at the decoupling , effects on structure formation are expected to be mild. 
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The plan of our paper is then the following. In Sec. II we will recall the basic equations and definitions with 
particular attention to the Alfven theorem and to the magnetic helicity theorem. In Sec. Ill we will introduce the 
class of configurations with finite energy and helicity. In Sec. IV we will introduce the gravitational effects of the 
knots and in Sec. V we will focus our attention on the problem of magnetic spikes in the gravitational wave spectrum. 
Section VI contains our concluding remarks. We made the choice of including in the Appendix different technical 
results which might be of help for the interested reader. 


II. DYNAMICAL SYMMETRIES IN HOT PLASMAS 


In a conformally flat metric of Friedman-Robertson-Walker (FRW) type the line element is 

= a^{r]) [di]'^ — dx '^], a{r]) ~ r]. 


( 2 . 1 ) 


Since we are going to exploit the Weyl invariance of the Maxwell fields in this type of backgrounds we stress, for our 
future convenience, that ( 77 , x) are the conformal coordinates which are related to the flat coordinates (t, y) by the 
relations dt] = dt/a{t), dx = dy/a{t). In the following we will denote by a prime the derivation with respect to 
conformal time and, when needed (for instance in Sec. Ill), by the over-dot the derivation with respect to the cosmic 
time t. 

Using the magnetohydrodynamical approximation the evolution equations of Maxwell helds at finite conductivity 
and diffusivity have to be supplemented by the Navier-Stokes equation. By exploiting conformal invariance they can 
be written as 


p + p ]v 


+ v- 'V 


p + p ]v 


+ vV 


p + p ]v 


= —Vp + J X H + 1 / 


-b ^V(V • v) 


H' + V X E = 0, V • i? = 0, 
yxH = J + E', V-H = 0, 

J = a^E + V X ) 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 


where v = x '{p) is the bulk velocity of the plas ma, a is the conductivity and ix is the shear viscosity coefficient; the 
curved space quantities appearing in Eqs. (^)-(|^) are related to the flat space ones by the metric rescaling 

E = a^S, H = a?B, a = acjc, p = Pc, P = of^Pc, J = a^j, v = aP'vc- (2-6) 


Notice that in Eqs. (|]^) V • A = 0 is not imposed as an additional equation on the fluid and indeed is not a valid 
equation for compressional waves which emerge from the theory |20 However, if the plasma is globally neutral 
on scales much larger than the Debye radius (i.e. Lo{T) ~ \JTj where rie is the average electron density ) the 
mean electric fields are effectively divergence free. After the epoch of e^-e” annihilation (i.e. T < 0.1 MeV) we have 
that Tie ^ 6.43 X 10“® Xe QbIt^qq T^, and it turns out that Li){T) ^ 1.9 x 10“'^(0.1 MeV/T) cm {xe is the ionization 
fraction and Qb the baryon density in critical units). If we compare the Debye scale with the horizon distance at 
T ~ 0.1 MeV we have that Lnite+e-) ~ 10^^ cm LD{te+e-)- Therefore, for our purposes the plasma description 
is indeed appropriate. 

The set of MHD equations can be studied with different “closures”. Some usual closures are the incompressibility 


(i.e. V • F = 0) and the adiabaticity, but other closures can be invented depending upon the problem at hand (like 
the isothermal closure or the closure J = constant). Since the evolution of the Universe can be described (within 
the hot Big-Bang model) by an adiabatic expansion, we can certainly require that p = Pc constant and that 
p = p/3. Moreover one can also assume that the fluid is incompressible and in this case, not only the magnetic 
lines will be divergence-less but also the velocity lines will share the same property. We notice, incidentally, that 
in the incompressible approximation we can also define a hydrodynamical helicity p7[ which allows the topological 
treatment also for the hydrodynamical part of the system. At the end of Sec. IV we will specifically consider a case 
where the incompressible closure is relaxed. 

In referring the reader to Appendix A for further details concerning the MHD evolution we want only to note that 
at high conductivity the Alfven theorem holds 


d 

dp 


H-dE = 


V X V X • FE, 
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(2.7) 















where S is an arbitrary closed surface which moves together with the plasma. If we are in the inertial regime (i.e. 
L > where is the magnetic diffusivity length), we can say that the expression appearing at the right hand side 
is subleading and the magnetic flux lines evolve glued to the plasma element. Moreover, the magnetic helicity TLm is 
conserved in the inertial range: 


'Hm = / cPxA 

Jv 




:f«" = -- 

dr] a . 


d^x H -V X H. 


( 2 . 8 ) 


Concerning Eqs. (|2.8| ) two comments are in order. In Eq. (|]^) the vector potential appears and, therefore it might 
seem that the expression is not gauge invariant. This is not the case. In fact A ■ H is not gauge invariant but, 
nonetheless, Hm is gauge-invariant since we will define the integration volume in such a way that the magnetic field 
H is parallel to the surface which bounds V and which we will call dV. Calling n is the unit vector normal to dV, the 
integral defining the magnetic helicity is gauge invariant provided H ■ n = 0 in dV. As we will consider finite energy 
fields in our examples of Section III we will have that H = 0 in any part of dV (which might be a closed surface at 


infinity). In Eq. (2.8) the term appearing under integration at the right hand side (i.e. i? • V x H) is sometimes 
called also magnetic helicity (or magnetic gyrotropy) and it is a g auge invariant measure of the diffusion rate of Hm 
at finite conductivity. Th e two theo rems given in Eqs. (|3)-(^ are reviewed and proven in Appendix A. 

Notice that Eqs. (^) and {2^) are truly dynamical symmetries of the plasma, in the sense that they are not 
inherent to static MHD configurations. Before ending this section we want to mention a further scale which turns out 
to be important in the discussion of the large scale effects of primordial magnetic fields, namely the magnetic Jeans 
scale B0. At t = tdec we have from Eq. ( |A.I2D that L^{tdec) — 4.7 x 10^° cm whereas the magnetic Jeans scale 
is, at the same epoch. 


\B{Uec)\ Mp 
Pc{i^dec) 


(2.9) 


In a matter-dominated Universe the magnetic Jeans scale evolves like the scale factor. In order to make our estimate 
more constrained we can assume a field sufficiently strong in order to rotate the polarisation plane of the CMBR at 
the decoupling [|| (i.e. \H{tdec)\ > 10“^ ) and we see that 

Lsjitdec) ^ 10^^ cm. (2.10) 


From this last numerical estimate we see that if the magnetic field is sufficiently strong, the magnetic Jeans length 
is of the same order of the horizon size (i.e. Lnitdec) LBj(tdec))- This result has interesting implications for the 
possible impact of magnetic field on structure formation. In particular, we have that if Lsjitdec) ~ Lniidec), then 
the density contrast induced by the Lorentz force term of Eqs. (2.2) and (A.17) can be of order 1 at the time of galaxy 
formation 0- In the opposite case (i.e. Lsjitdec) ^ Lniidec)), from the Lorentz force term we cannot expect any 
significant effect on structure formation. In Sec. V we will come back to this problem. 


III. MAGNETIC KNOT CONFIGURATIONS 


Since both magnetic helicity and magnetic flux are conserved during the dynamical evolution in the inertial range, 
configurations with non trivial topological structure present at some initial time will be preserved all along the 
dynamical evolution. 

In this Section we will provide some examples of magnetic knot configurations. In order to perform various calcu¬ 
lations, it is useful to employ configurations which are non singular at large and small distances. In this way all the 
integrals defining the helicity, the gyrotropy and the total energy, will be automatically well defined. This is exactly 
what happens in the following examples. Consider the magnetic field with (spherical) components 
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FIG. 1. We plot the moduli of the poloidal (right) and toroidal (left) components of the held conhgurations reported in Eq. 
([l.l|). In these topographic maps dark regions correspond to low intensity spots of the held whereas brighter regions correspond 
to higher intensities. The toroidal component does not depend on n. The poloidal component is instead computed in the case 
of n = 100. Moreover in both case we took \Ho{tdec)\ ~ 0.001 Gauss. Notice that in the limit of n —> 0 the poloidal component 
goes to zero. For r = L/Ls >2 the held intensity is suppressed. 



In Eq. represent the typical scale of the seed and TZ = r/Lg. Moreover Bg = HqL^ is the dimensionless 

amplitude of the magnetic field; n is just an integer number. It can be also intuitively useful to report explicitly the 
moduli of the poloidal and toroidal components of H 


\Hp{n, 0 ,n)\ 


-T- 




8Bo TZsin0 


(3.2) 


where the poloidal and toroidal components of the field are defined in the standard way {Hp = HrCr + Hgeg, 
Ht = H(j,e^). The configurations ( fl.lD were firstly introduced in in the context of the study of the topological 
properties of the electromagnetic flux lines. 



whereas its length is proportional to the held intensity. We notice that for large n the held is concentrated, in practice, in the 
core of the knot (i.e. r/La < 1). This can be also argued from Fig. (|^. 


From the explicit expressions of the vector potentials 
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the magnetic helicity can be computed 


Hm 


A ■ Hd^x = 



SuBq TZ^dTZ 
[7^2 + l]‘‘ 


= nBl 


(3.4) 


Notice that the integration is perfectly convergent over the whole space and, since the field is i? = 0 in any part of 
dV at infinity, this quantity is also gauge invariant. We can compute also the total energy of the field, namely 



(3.5) 


The magnetic configurations we just presented are solutions of the MHD equations (in the inertial range) with large 
Prandtl number and small [3 parameter (/3 «C 1, Ptm ^ 1) for a velocity field v = ±H/[p + p\ (see Appendix A 
for more detailed explanation of these notations). 

In the limit n —*■ 0 the total helicity goes to zero and, as a consequence, the magnetic field is purely toroidal. In the 
case where the helicity is non vanishing we can also compute the volume integral of the magnetic gy rotr opy H - V x H 
which, according to Eq. (^), measures the diffusion rate of Hm at finite conductivity. From Eq. (|3.l|) we have that 


H -V X Hd'^x = 


256 Bl n r°° n^dn 
Jo (T+W 


57Ym 


(3.6) 


Using this last result into Eq. (|2.8| ) we obtain the explicit version of the helicity dilution equation valid in the context 
of our configurations, namely 




(3.7) 


IV. GRAVITATIONAL WAVES AND MAGNETIC KNOTS 

Magnetic knot configurations present right after the epoch of electron-positron annihilation can act as seed sources 
for the evolution equations of the scalar, vector and tensor metric fluctuations. In fact, the spatial part of the knot 
energy momentum tensor (computed in Eqs. ( |B.5| ) of Appendix B) is a symmetric three-dimensional tensor of rank 
two with six degrees of freedom which transform differently under three-dimensional rotations on the 77 = constant 
hypersurface. More precisely the spatial components of the knot energy-momentum tensor can be decomposed in 
a trace part plus a traceless scalar with one (scalar) degree of freedom each, and vector and tensor parts with two 
degrees of freedom each. The energy momentum tensor of the knot has no homogeneous background contribution. 
This fact has two important consequences. On one hand the knot does not affect the background evolution of a 
radiation (or matter) dominated Universe, on the other hand the energy momentum tensor is automatically invariant 
for infinitesimal (gauge) coordinate transformations. In the following discussion we will focus our attention on the 
evolution of the tensor modes. The tensor modes of the metric correspond to gravitational waves . 

Consider then a perturbation of a (homogeneous and isotropic) conformally flat metric 

^ 5m--(??)+% i.(u??)- (4.1) 

whose fluctuating part, (5< 7 ^ 1 /, contains scalar, vector and tensor modes 

= hg(,'^^(f,77) -k(55^^^(f,77) -k5g^'^^(f,77), (4.2) 

The corresponding equations of motion will have, as source term, the scalar vector and tensor modes of the knot 

energy-momentum tensor 

STfj,^{x,r]) = STl^f,\x,ij) + 5Tli^\x,r]) + 5Tli^\x,T]). (4.3) 

Tensor perturbations of the metric are constructed using a symmetric three-tensor which satisfies the constraints 

hi = = 0, (4.4) 
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which implies that hij does not contain parts transforming as scalars or vectors. Therefore the perturbed tensor 
components of the metric will be 

hlo = 0, (5^0^ = 0, Sgij = -a^(r])h,j{x,ri), (4.5) 

and the perturbed line element becomes 






drj — {'jij + hij)dx^dx^ 


(4.6) 


where 7 ^ denotes the spatial part of the background metric. The evolution equations for the tensor modes of the 
geometry can then be written as 




(4.7) 


('T'\ T 

wher e 5G\J are the tensor components of the Einstein tensor G^jy = — ^gfj,uR perturbed according to Eqs. 

(4^)-(^). Eq. (^) equation becomes, in conformal time, 


h" + 2Hh[j - = -WttGtIP. 

Notice that the energy-momentum tensor can be decomposed, in Fourier space, as 
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3^' 


e Ty(A;)d'^A:, T(fc) = ni(k), 
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+ ( hkj - ) + ( kiTp\k) + kjTp\k) ) -b ^(fc), ki = 
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where 
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(4.9) 


T^^Hk) = -j^khT^Ak) - 

= Tim{k)km - kkjTjm{k)km, 

Tpkk) = Tijik + ^ ( kkj - Stj]T{k) + ^ ( kikj + Sij ) kmknTmn{k) - kikmTmj{k) - kmkjTim{k) 


(4.10) 


In Appendix C the Fourier transforms of each component of the energy momentum tensor is reported. Here we 
will simply discuss the results of the expansions in the two limits relevant for physical applications. For frequencies 
much larger than the typical frequency of the knot (i.e. fc > fcg ^ have that the Fourier transform of each 

component of the energy momentum tensor is exponentially suppressed exp[—fc/fcg], whereas for small frequencies 
(infra-red limit) the various components have different power-law behaviours: 
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(4.11) 


( O denotes the Landau symbol; notice that, within the square brackets we gave the infra-red expansion of the Fourier 
transform, whereas, outside the brackets we kept the leading [ultra-violet] e xpon ential suppression which always 
factorizes according to the exact results reported in Appendix C). From Eqs. ( [l.ll| ) we can see that the off-diagonal 
terms are alw ays subleading (at large scales) if compared with the diagonal ones. 

From Eq. ( |4.8| ) it is possible to define the (Fourier space) evolution of the two transverse and traceless degrees 
(metric) of freedom which reads 


hij (f, 7 ) = 

Tp\x,v) = 


(2^)t 


(27r)t 


<Ak qAhQ{k,r])e^'^'^ + qlj’hjj{k,r])e^'^'^ 


(4.12) 
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with 




ef\k)e^;\k)-ef\k)ef\k) 




e\^\k)ef\k) + ef{k)e^^\k) 


(4.13) 


Notice that k/\k\, e(i), 6 ( 2 ) are three unit vectors orthogonal to each others: 
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\k\ 

6 ( 1 ) = 


= ( kx,ky,kzj = ^sin 0 cossin 0 sincos 6 * 

kx ky 


1-2 I h 2 ^ h '2 4 - h '2 

'^x^ '^y Y w ftj, 


, 0 = sin((), — cos (j),0 ], 


6 ( 2 ) = 


kx kz 


1 + kl \ = ±( cos 9 cos cj), cos 0 sin ((i, — sin0 ), 


t-2 _1_ 1.2 ^ / 62 I r,2 


(4.14) 


(the ± in the last formula refers to the cases when 9 > 7 r /2 [+] and 9 < 7 r /2 [-]). After having extracted the tensor 
modes from the Fourier space components of the energy-momentum tensor we can find, with some algebra, T 0 (fc) and 
r^ik): 


TT Bn 


Te(fc) = J 2 + K) 2(fc): r^ik) = -J 2 r- + + ‘^^z) Q{k) 


TT Bl 


Q{k) = 


7n^ - 5 
40 




2 \ 1 


(4.15) 


Then the evolution equations for the two polarizations become 

/i0 -I- 27fft.0 -I- k^hf^ = — 167rG'T0, 
h'^ -I- 2Hh'^ + k^ht^ = -IbTrGr®. 

We solve the evolution equations for /i 0 and ft.® by requiring that at rji 

h®{m) = *©(fti) = 0, = o, 

with the result that 


(4.16) 

(4.17) 


h®{k,ri) = ^^^r©(fc) 
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krii cos [fc (?7 — Yi)] + sin [k{r] — yi)] — krj 
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(4.18) 


Using these expression we can compute the energy density of the produced gravitational waves. Up to now we 
were dealing only with one configuration. In principle after the epoch of electron-positron annihilation different 
configurations can be present. If the configurations are not correlated the polarisations of the produced gravitational 
waves will also be stochastically distributed, which means that 


{<ltf ik)q"(i')ik')) = Sss'S^^\k-k'), 


ilj V”vr(s')V"' 7/ — — (4-19) 

where < ... > denotes a stochastic average. The total energy density radiated in gravitational waves will then be 
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Define now the comoving frequency to = k/a and recall that krj = lo/H (where H = d/a is the Hubble factor in 
cosmic time). By taking rji around the epoch of electron positron annihilation (i.e. T ^ 0.1 MeV) we can compute 
J-{kr]) at any interesting time rjQ as a function of the frequency. Taking rjQ ~ r]ciec (when T ~ 0.26 eV) the (present) 
decoupling frequency will be tUdec ~ lO^^s whereas the frequency corresponding to the present horizon will be 
loq ~ 3.2 X The present frequency corresponding to T ~ 0.1 Mev is ujg ~ 10“4i Hz and for the range 

LOdec < W < Ws we have that \T{kridec)\ = \^{‘^l^dec)\ ~ Y^l^dec)~‘^- 


V. MAGNETIC SPIKES IN THE GRAVITATIONAL WAVE SPECTRUM 

For the discussion of the gravitational wave spectra produced by magnetic knots, it is useful to define the energy 
spectrum in critical units, 


Dgw(w, t) = 


1 dpGW 
Pent d log W ’ 


(5.1) 


which will allow a comparison with t he spectra produ ced, v ia gr avitational instability, in the context of ordinary 
inflationary model. Using Eqs. ( 4.15 ) into Eqs. ( 4.20| ) and ( 4.21 ) we can perform the angular integration over the 
momenta and we are left with the integration over dlogfc. Passing to the physical frequencies and using Eq. (|5.1| 
get that 


we 
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OJdec < CO <UJs 


(5.2) 


(notice that A = Hq/p = BqY)/ pc{t) and D.y(t) is simply the fraction of critical energy density present in the form 
of radiation at a given observation time t). For the purposes of the present Section it is also useful to introduce the 
logarithmic energy spectrum for the knot energy density 


t) 


1 dpKN 
Pent d log W ’ 


(5.3) 


where Pkn(Y) is given by 'roo(fc) after the integration over the directions of k. 
that, in the range ujdec < uj < Us, 


Using Eqs. (C.l) and (C.2) we get 


+ 1)A 


3 r 


1 + 0 


UJ 


(5.4) 


Obviously, for the consistency of our considerations we have to impose, for all the frequencies and for any time, that 
OiCAr(a;,t) < 1. This simply means that the energy density of the knot configurations does not modify the (radiation 
domi nate d) background evolution and can be always treated as a small perturbation. Given the analytic form of 
Eq. (5.4) we see that the largest contribution to Q,KN{+,t) comes from modes of the order of the size of the knot. 
For higher frequencies the exponential damping becomes effective. Recalling that, today, — 10“^, the critical 

energy condition Q,KN{+>,t) < 1 implies that, for ui < lOs, 


logio ^ + logio(»^^ + 1) < 4.53. 


(5.5) 


It is also worth mentioning that in the above region of parameters VLQw{^,t) < 1 is automatically satisfied. 

We want now to compare the logarithmic energy spectrum of the gravitational waves produced by magnetic knots 
with some (possible) inflationary spectra. It is indeed well known that any transition in the curvature inevitably 
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amplifies the tensor fluctuations of the metric ||30|-|32|. In particular the transition from a primordial inflatio nary 
phase to a decelerated one is associated with the production of a stochastic background of gravitational waves [Q, 
which decreases as for loq < oj < ujdec and stays flat for Udec < w < Uds ( ^ds = Hz is the 

present frequency associated with the maximal curvature scale Hds reached during inflation). The gravitational 
waves spectrum induced by the inflationary transition can then be written as 




to) = ^-f{to) 


Hds 

Mp 

Hds 

Mp 


O^dec ^ UJ < OJdS^ 


^dec 


LJO <UJ < UJdec- 


(5.6) 


There are various bounds which should be taken into account. At large scale the most significant one comes from the 
anisotropy of the CMBR which implies that 


^Gwi^,to)hfQQ < 7 X 10 , for w ~ loq. 


For uJdec < Ld < ujds, using the bound (^) into Eqs. ( |5.6| ), we have 

to) < 10 “^^, 

or, in terms of ffdS P^ , 

HdS 


-6 


At intermediate frequencies u 


Mp 

10“® Hz, the extreme regularity of the pulsar’s pulses 
VLcw{<-op,to) < 10“®, for uip ~ 10“®. 


imposes 


(5.7) 


(5.8) 


(5.9) 


(5.10) 


Due to the flatness of the logarithmic energy spectrum of Eq. ( |5.6D around wp, we can clearly see that the pulsar’s 
bound is easily satisfied if Eq. (5.8) is satisfied. 


We want now to compare the spectrum given in Eq. (5.2) with the one produced thanks to the inflationary 
(param etric ) amplification of the quantum mechanical fluctuations of the geometry reported in Eq. ( |5.6| ). If we look 
at Eq. (^) around ujg we can see that the gravitational waves produced by the tensor modes of the knot can exceed 
10“^^ (i.e. the maximal allowed amplitude of the inflationary logarithmic energy spectrum in critical units) provided: 


logio A + - log 


10 


/7n2 - 5 
20“ 


> -0.537. 


(5.11) 


The last equation was simply obtained by requiring that Q,Gw(oj,t) given in Eq. (5.2) can be larger than 10“^^ for 
L Odec < dj < ids- In Order to be consistent with our approximations we have, therefore to impose simultaneously Eqs. 
( ^.11 ) and ( ^.5[ ). The results are reported in Eig. || in terms of the two parameter of the model A (the intensity of 
the magnetic field in critical units) and n the number of knots of the configuration. 



parameters A and n necessary in order to give a signal of the same order of the stochastic backgrounds of inflationary origin 
(lower line). For values of the parameters within the shaded region, the gravitational wave signal produced by the tensor 
part of the knot energy-momentum tensor exceeds the inflationary prediction also of four orders of magnitude in the (present) 
frequency range ids ~ Hz. More quantitative illustration is given in Fig. W. 
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In Fig. ^ there are, formally, two parameters. One is A and the other is n. From the physical point of view, 
A represents the energy scale of the configuration in the limit n ^ 0, whereas, n is exactly the magnetic helicity 
discussed in the previous Sections. We stress the fact that the helicity is conserved for a wide interval of scales after 
tg+e-, and, therefore, it represents a very good classical “label” which can be used in order to find the gravitational 
“imprints” of magnetic knots. 



(co/Hz) 

FIG. 4. We report various spikes in the gravitational wave spectrum computed for different values of n. The values of A 
is 3 for the two lower spikes and it is 1 for the upper one. In this plot we reported the expected signal computed assuming 
that quantum mechanical fluctuations were amplified by the transition between the inflationary phase and radiation dominated 
phase. Notice that the exponential damping ensures the compatibility with the pulsar bound. 


The spectrum of gravitational waves induced by the knot of Eq. (5.2) decreases as (lo/lus) for to < tos and it is 
exponentially suppressed for w > Wg. This peculiar behaviour has two important consequences. Firstly, the exponential 
suppression guarantees that the pulsar bound (see Eq. ( 5.101 )) does not constrain the scenario since at that frequency 
the gravitational signal of the knot will be suppressed of a factor of the order of exp [—uip/ujs] ~ exp (—10^). Secondly, 
the steep infrared behaviour guarantees a quick suppression at scales even larger corresponding to frequencies of 
the order of 10“^® Hz. Therefore, the signature of these configurations would be compatible with a spike in the 
gravitational wave spectrum for frequencies of the order of Wg. Various magnetic spikes are depicted in Eig. We 
can clearly see that a large magnetic helicity certainly helps in getting a signal which can be substantially larger (four 
orders of magnitude in View) than the one provided by the stochastic backgrounds of inflationary origin in the same 
interval of frequencies. In Fig. ^ we report also the sensitivity of the Laser Interferometer Space Antenna (LISA) for 
typical (present) frequencies of the order of ujl ~ 10“^ — 10® Hz We would also like to notice that the spikes 
we just described are also well below the, so-called, nucleosynthesis bound |^. In fact the simplest (homogeneous, 
isotropic) big-bang nucleosynthesis scenario imposes and (indirect) constraint on the gravitational wave spectrum 
which can be written as 


J d\ogujVlGw{i^,to) < 0.2 Vlj{to). 


(5.12) 


This constraint implies that VtGw{^,to) should be below 10“® as it occurs for all the modes of our spectrum. 

Magnetic knots present right after the epoch of electron-positron annihilation might have, in principle, some in¬ 
fluence on the process of galaxy formation. Our discussion will be quite qualitative and will essentially follow the 
approach of Ref. [^. Consider the model of a pressure-less fluid made of baryons and of dark matter particles 
with unperturbed velocity fields Ug = (l/a,0), = (l/a,0). The velocity fluid perturbations will then be given, 

generically by Su^ = n*/a. Since in our discussion we assumed from the very beginning that (3 < 1 (as implied by 
the critical energy condition on the energy density of the knot) the pressure gradients will also be subleading and, 
therefore the system for the perturbed density inhomogeneities will be (in the absence of viscous terms) 


6p -|- Hd'g — 


s'd + hSd - 


VIrSb + VIdSd 
VIbSb + VId5d 


V • X V X iJ] 
[pd + Pb] 

= 0 . 


(5.13) 

(5.14) 
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(5.15) 


where Ti = a'/a is the Hub ble f actor in conformal time, 6b = 5b/pb, 6d = 6pd/pd and VIb = Pb/[pb + pu], 
= Pd/[pb + Pd]- In eq- (5.13) the contribution of the Lorentz force has been written by taking into account the 
fact the electric field are suppressed at high conductivity. Multiplying Eq. ( |3.13| ) by Hb and Eq. ( 5.14 ) by Hb we 
find, after having summed up the two obtained equations, that the total density contrast 6 = flB6B + ^d6d obeys a 
decoupled equation 


„ , 3 n VJHxVxHI 

5" + m' - -n^s = Qb -r- 

2 [pd + Pb\ 


Notice that, in cosmic time this equation becomes 


5 + 

a 2 


S — Qt 


V • [H X V X 
[o2(pb+Pb)] 


(5.16) 


(5.17) 


The solution of this equation in a matter dominated Universe after recombination (t > tree) can be obtained imposing 
the initial conditions 6{tree) = 0 and 6{tree) = 0 with the result that at a given length scale L the density contrast is 
given by 


6{L,t) 


144 

^Bj (^rec) 

2 

a{t) 

157r3 

Lsij'rec 


Ctitrec') _ 


Q{\n) 


(5.18) 


where A = L/Ls and Q{X,n) contains the spatial dependence. In the case of the configuration given by Eqs. 0) 
0(A, n) turns out to be: 


0(A,n) =Po(A) + P„(A), 


Po{X) = 
Pu{X) = 


4 - IQA^ + 19A^ + (27A2 - 15A^) cos26» 


2A® - 16X^ + 19A2 - 5 - am 28 - 13A2 + A^ cos 29 


(5.19) 


Notice that Q{X,n) receives its maximal contribution around A = 1. This simply means that the effect of the Lorentz 
force term becomes important only for scales comparable with the size of the knot which, in our scenario turns out 
to be of the order of the horizon size right after the epoch of electron-positron annihilation. 

In order to affect structure formation through the Lorentz force we need that, today 


pBj{to) ~ 0.03 Ls{to) 


(5.20) 


with Ls{to) ~ 1 Mpc. From this last condition follows that H{tdee) 
decoupling time we have that 


10 ^Gauss. Blue-shifting Eq. 


( 5.20| ) at the 


pBjitdee) ^ 2.7 X 10^® cm. 


(5.21) 


Now, the typical scale of the knot is Ls{te+e-) ~ Lnite+e-)- This means that Ls{tdee) ~ 10^®cm. Thus we have that, 
in our scenario, Lbj ~ 10 Ls{tdee)- As we just said, in order to have a significant impact on structure formation, 
we would need Ls{tdee) ^ Ls{tdee), whereas, in our case, we have Ls(tdec) 0.1 Ls(tdec) which might have some 
mild impact on structure formation. In the spirit of the order of magnitude estimate presented in the last part of the 
present Section we notice that, if appreciable effects occur on structure formation (i.e. \H{tdee)\ ~ 10“^Gauss) the 
very same amplitude of the magnetic field can depolarise the CMBR . 


VI. CONCLUSIONS 

The dynamical symmetries of hot plasmas impose the conservation of the magnetic flux and of the magnetic helicity. 
In this paper we gave some examples of magnetohydrodynamical configurations carrying a finite magnetic helicity. 
We argued that, thanks to the helicity conservation the topological properties of the divergence free magnetic field 
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lines are the only dynamically motivated parameter which could be used in order to discuss the possible implications 
of magnetic fields in the early Universe. 

Based on the fact that after the epoch of electron-positron annihilation the MHD are fully valid we argued that if 
magnetic knots were present they could have acted as a source for the evolution equation of the scalar, vector and 
tensor modes of the metric inhomogeneities. We focused our attention on the tensor modes and we found that the 
induced gravitational wave spectrum is sharply peaked around a (present) frequency of the order of 10“^^ -10“^^ Hz. 
The typical amplitude of the spike depends upon the helicity of the configuration or, more physically, upon the the 
number of twists and knots in the magnetic flux lines. The amplitude of the spikes can be substantially larger than 
the gravitational waves signal provided, in the same interval of frequencies, by inflationary models. 

Many interesting points are left for future discussions. First of all there are technical points which could be 
addressed. In this paper we only took into account the tensor modes. Of course it is certainly possible to extend 
the above discussion to scalar and vector modes. It is not excluded that some interesting consequences for the small 
angular scales anisotropies can be addressed. Moreover we would like to stress that the present results hold in the case 
of one configuration. It would be interesting to generalise these results to more complicated statistical distributions 
of magnetic knots. This problem has not been addressed in our present investigation. 

There are also deeper questions. In this investigation we were dealingwith ordinary MHD. It is, in principle possible 
to extend our discussion to earlier times (like the electroweak epoch) ill- 
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APPENDIX A: MAGNETIC FLUX AND HELICITY THEOREMS 


Consider and arbitrary closed surface E which moves with the plasma. Then, by definition of the bulk velocity of 
the plasma {v we can also write dS = v x dl drj. The (total) time derivative of the flux can therefore be expressed as 


■^ [ H ■df; = [ ^ ■d^+ [ Hxv-dl 

dr] Jy, ov JdY. 


(A.l) 

where 9E is the boundary of E. Using now the properties of the vector products (i.e. H x v ■ dl = —v x H ■ dl) we 


can express v x H though the Ohm law given in Eq. (2.5) and we obtain that 


V X H = —E - 


1 , 


-V X H 


Using now Eq. (AU) together with the Stokes theorem we obtain 




dH 

dr] 


+ V X U 


•dE- - f VxVxid-dS 


(A.2) 


(A.3) 


From Eq. (2.3) given in Sec. II the first part at the right hand side of Eq. (A.3) is zero and Eq. (2.7), expressing the 
Alfven theorem, is recovered. 

With a similar technique we can show that the conservation of the magnetic helicity also holds. Consider a closed 
volume in the plasma, then we can write that dU = d^x = v± ■ dE dr] = n-v± dE dr] where n is the unit vector normal 
to E (the boundary of V, i.e. E = dV) and Tl is the component of the bulk velocity orthogonal to dV. The (total) 
time derivative of the magnetic helicity can now be written as 


■^Hm= [ d^x^[A-H] 
dr] Jv or]^ 


A ■ Hv± ■ fidTi 


ldV=Y, 


(A.4) 


By now writing down explicitly the partial derivative at the right hand side of the previous equation we can use the 
MHD equations (recall that E = — A' and that H = W x A). Finally using again the Ohm law a nd tr ansforming the 
obtained surface integrals into volume integrals (through the divergence theorem) we get to Eq. (2.8). 

We want to stress the fact the the magnetic helicity is indeed a gauge invariant quantity under the conditions stated 
in Sec. II. Consider a gauge transformation 


A ^ A + 


(A.5) 


then the magnetic helicity changes as 


[ d^xA-H^ f d^xA-H+ [ d^x^-[^H] 
Jv Jv Jv 


(A.6) 


(in the second term at the right hand side we used the fact that the magnetic field is divergence free). By now using 
the divergence theorem we can express the volume integral as 


d'^xV 


[fH] = / 

Jd 


■ fidTi 


(A.7) 




Now if, as we required, id • n = 0 in dV, the integral is exactly zero and TLm is gauge invariant. Since in Sec. II 
we considered fields going to zero at infinity id = 0 in any part of dV at infinity. Notice that the field go to zero 
at infinity is not really necessary. In fact we could just define the integration volume by slicing U in a collection of 
closed flux tubes. Now by definition of flux tube id • ri = 0 if n is the unit vector orthogonal to the flux tube. 

In order to give a self-contained presentation of our problem we want to recall the main limits where the MHD 
equations can be discussed. This part of the Appendix represents an ideal follow up of the considerations presented 
in Section II. 


Using the Ohm law of Eq. (2.5) in order to express the induced electric field in terms of the current density and 


taking into account that J = V x dd, we obtain from Eq. (2.4) the standard form of the magnetic diffusivity equation 


Id' = V X [vx H]+ --V^dd, V • id = 0. 


(A.8) 
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Assuming the incompressible closure together with the adiabatic assumption (i.e. p' = Q, p' = 0) we get a simpler 
form of the Navier Stokes equation 


'' + [u- Vjt) = — 


1 


P + P 


\H[ 


P' 




P + P 


-V^v. 


(A.9) 


In writing this equation we used the fact that 2 J x H = —V|i?p + 2[H ■ V]i?. Eqs. ( |A.8| ) and (AJ) contain the 
magnetic and thermal diffusivity scales. We will focus our attention on the case of an electromagnetic plasma whose 
temperature is T ^ 100 GeV (in the opposite limit [T > 100 GeV] the equations appearing in the present Section 
should indeed be generalized to the case of finite chemical potential leading to the anomalous MHD equations |^). 
An estimate of the plasma conductivity is necessary for the calculation of the magnetic diffusivity scale. Before the 
epoch of e+-e“ annihilation (i.e. T > 2me ~ 1 Me V) th e conductivity can be written as Uc = jci&ra where A is a 
constant of order unity (recall also that, from Eq. ( |2.6| ), a{r]) = CTc a(ry)). For temperatures smaller than the weak 
interactions equilibration temperature (i. e. T < 0.2 irie) the electron velocity can be roughly given as v = \/T/rrie 
whereas the collision cross-section is Ocoii — {demlT)'^ (up to the Coulomb logarithm which is of order one in our 
case). Therefore the conductivity will be simply estimated by 


(A.IO) 





From Eq. (A.S) the magnetic diffusivity scale will be 

Loiv) 


V 


aca{r]) 

Using Eq. ( A.lOj ) we can say that, after the weak decoupling, 

T 


L„{T) ~ 4.3 X 10 


10 


0.308eV j 


cm, T < 0.2 rUe 


(A.ll) 


(A.12) 


26| | range of the magnetic field spectrum where the effect of the finite value of the conductivity 
y neglected (the so called ideal limit of the MHD equations). The thermal diffusivity coefficient 


(where we used as reference the recombination temperature T^ec 0.308 eV). Notice that this scale is incredibly 
small if compared with the horizon size at each corresponding epoch. The magnetic diffusivity scale tells us which 
is the “inertial” |^, 
can be approximate 
is given by v/ (p -b p) ~ 0.2 X-y{T) where 

X^(T) - 9.106 X 10^° x-^ (fls hlQo)~\XleY/TfT-^ (A.13) 

( Xe is the ionization fraction and XIb hfpQ is the fraction of critical energy density of baryons). Therefore the thermal 
diffusivity scale will be 


Lth{T) ~ 1.5 X IQ^^Xe ^ (Dsh^oo)”^ f ' T-\ 


(A.14) 


It is important to notice that the ratio of these two quantities gives essentially the (magnetic) Prandtl number 


Pr^{T) = 3.5 X lO^^ x-\VLBhl^^)-^ '. 


(A.15) 


Since Prm{Tdec) ^ 1 the effect of the Lorentz force terms in Eq. (A.9) cannot be neglected. Another interesting 
dimensionless ratio is what is usually called /3-parameter ||^ 


m = 




45 \H{T)\'^ 


2p 2TT^geff 


(A.16) 


If /3 <C 1, then the term V[p + |i7p/2] can be neglected. The critical energy condition applied on the magnetic field 
backgrounds at any given temperature in the radiation dominated epoch enforces by itself this assumption. In order 
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to give a more quantitative estimate, let us suppose that the magnetic field is strong enough to produce an appreciable 
depolarisation of the CMBR. If this is the case we should have \H{tdec)\ ^ 10 ~^ G a uss around Tdec ^ 0.26 eV. 
This means that (3{Tdec) ~ 2.7 10“^. Therefore for large Prm and small (3 Eqs. (A.8)-(A.9) can be further simplified: 


1 


77' = V X [vx H]. 


(A.17) 

(A.18) 


APPENDIX B: ENERGY-MOMENTUM TENSOR FOR KNOT CONFIGURATIONS 


In this appendix we want to give the components of the energy momentum tensor which are used in our calcula¬ 
tions. We want to express the energy-momentum tensor in Cartesian coordinates where some estimates become more 
tractable. Let us start by giving the cartesian components of the magnetic field: 


H4X) = 

Hy{X) = 

H,{X) = 


4Bo 2y-2 n XZ 

[1 + ^-2-7 3^2+ 2 : 2 ]^’ 
4Bo 2X + 2 nyZ 

[i + A’2-t3;2 + 22]^’ 
4Bo n[X‘^ +y‘^ -Z^-l] 

[i + ^2_^3;2_^22]^’ 


(B.l) 


where, as usual, X = xjhs. y = y/Lg, Z = zjLs. The spherical components reported in Eq. ( [l.l| ) are simply related 
to the Cartesian ones. In Cartesian coordinates we have 


H{X) — Hx^X + HyCy 




(B.2) 


where Cx, Cy and are three mutually orthogonal unit vectors. The three unit vectors in spherical polar coordinates 
are determined by computing the Jacobian matrix with the result that: 


Gr - Gnr 

n 


y ^ 

Gy 

n ^ 


n 


zx 


eg = 


G(j) — 


7^^/7^2 - Z2 
3^ 


Z3^ ^ v^^2Tr^^ 


7^\/7^2 - Z2 
A” 


y/X‘^ + 3^2 


(B.3) 


The polar expression for the components of the magnetic fields are obtained, from Eqs. (BT), (B^) and (B^) as 
Hr = H ■ Br, Hg = H ■ eg, = H ■ etj,. In order to get to recover the result of Eq. (^) the well known relations 
between spherical and Cartesian coordinates should be used (i.e. X = 72. sin 0 sin 0, y = 72. sin 0 cos </), Z = 72. cos 0). 

Recall now that H{X) = a?B. In terms of B the Maxwell field strength is given by Fij = o^eijkB^ and the 
corresponding energy-momentum tensor will be 


5Tfj,r, = ( - -gyuFafj 


paP 


(B.4) 


(we kept the notation STy_i, used in Section 4 in order to stress the fact that B is nothing but a fluctuation on the 
homogeneous background). It is often useful to express the components of the energy momentum tensor (as we did) in 
terms of the rescaled fields H. The relation between the energy-momentum tensors is = a~^Ty,i/ {ST^'^ = o“®r^'^), 

where now t^i, is only expressed in terms of H (remember that we work always in conformal time ). The components 
of the Energy-momentum tensor are: 

8B^ 4{x'^ + y^) + n'^[i + x"^ - 2y^ + y^ + 2Z^ + 2y^z^ + z'^ + 2x'^{y^ + z^ - 1)] 
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TxxiX) = 


‘yy 


{X) = - 


8B^ 4(X^ - y^) + 16 nxyz + n^[l + X^- 2y^ + 3^^ + 22 ^ + 2Z^y^ + Z^ + 2X‘^ {y"^ - 3Z^ - 1)] 

[1 + ^ 2+372 + 2 : 2 ]® 

8^2 4(A’2 - y^) + 16xyz -n'^[l + X^- 23^2 + 374 + 2 Z 2 - 63 ^ 2^2 + + 2 ^- 2(372 + Z 2 - l)] 


.XX) = - 


LW [1 + ^2 + 372 + 22 ] ® 

8^2 - 4 (A ’2 + 372) + n2 [1 + ^”4 - 23^2 + 374 + 2Z2 - 63^2^2 + Z 4 + 2^-2(372 - 3 Z 2 - 1)] 




_ 6iBUy-nXZ){X + nyZ) . 

Txy(X) — r, , ,n „on 6 ’ Txz\X) — 


[1 + ^-2 + 372 + Z 2 ] ® 

32Bln (-3^ + nXZ) {X^ +y^-Z^-l) 


TyziX) = 


LW [1 + <^2+372+Z2]® ’ LW 

32Bln {X + nyZ) {X^ +y'^-Z^-l) 


[1 + ^-2 + 372 + Z 2 ] 


Ll7r2 


[l + A’2 + 372 +Z2]' 


(B.5) 


Concerning this energy-momentum tensor few comments are in order. We did not include any electric field. The 
reason is simply that all the parts of the energy-momentum tensor containing electric parts are sub-leading at high 
conductivity. In fact, as we discussed in Section II in our approximations the electric fields are given by i? ~ V x i7/cr 
and, therefore, negligible in the ideal approximation of MHD. However, if one would like to study precisely the border 
region of the i deal approximation (i.e. typical scales of the knot comparable with the magnetic diffusivity scale given 
by Eq. ( A.12| )). then, the electric components should be included (sometimes this approximation scheme is named 
resistive MHD [Q). In our discussion the knot scale Lg is comparable with the magnetic Jeans length and then much 
larger than the magnetic diffusivity scale. There are interesting limits where the various components of could be 
investigated like n ^ 0. 


APPENDIX C: FOURIER TRANSFORMS OF THE ENERGY-MOMENTUM TENSOR COMPONENTS 


In this Appendix we report the Fourier transforms of the components of the energy momentum tensor of our 
configurations which are crucial for the discussion of the induced inhomogeneities. 




TT^Lfk^ 

+ —IXk) 

240 ^ 


-I- I4(k) 

80 ^ ^ 


2(n^ - 2 )kI -f 2{rX + 2 )kI - 2n^K^ -f 7n 


—riK^ — n{5Ky + k 1 + SiKxKyK^ 


n 

iMo 


hiK) 


+ 2nl{nl-3nX) + {nl + nl) 


y{k) = 




1 


hiK,) - Thin) 




240 


5‘ 

( 2^2 + 2^2 _ 5 k ^) — SlKxKyKz 




2[rX X 2)«;2 -|- 2{rX — ‘2)Ky — 2rXK\ -|- 7n^ 
2 


-h 


n 


-J2{k) 


, + 2s;2(k2 + ^2) -I- (^4 _ QKyKl + kX) 




5760 

2n^K^ — 2(n^ -I- 2)kx — 2(rX + 2)Ky 


+ ^/2(«) kI-2{kI + 3kI) -^M^) 


1920 

64i?g { . 2 1 ^ \ ^ ^ \ 




32^2 


M2o"4») 


1 


/ 2 2 I 2 \ 

KxKz(Hx ~ + "nKy) 


480 

i 


160 


/4(k) 


2nKxKz — 3iKy 

1 


- T7I^5(^) 


32 H 2 f . ^ 

^ ^2^4^! 1 1920 ^2 (k) 


.2 .2 , „.2 


/ 2 2 I 2 \ 

KyKz{Ky — Kz,+ riK^) 






‘Xny + — 3kz — Kynl + Kx{ 8 — iriKz) 
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Too[k) = 


1 


7r2Lf/cf [ 1920 
1 


/2(k) 


“t“ ^z') ~k ^x^y ~k ‘^^x^z 


+ 


4 — 77.^ 


h{K) + n"^ Ie{K) 


240 

1 


hiK) 


3n^(/v^ + '^y + ^z) + + 2^z) 


- 80 ^‘<''> 


1571 ^ - 2n^{Kl + kI- kI) - 4 ( k 2 + nl) 


(C.l) 


Notice that = /c^/fes and that k = \/k^. In the same way as we reported the Cartesian components of the energy 
momentum tensor in terms of the “rescaled” coordinates (zf, y, Z) = (x, y, z)/ Lg^ also dealing with the corresponding 
Fourier components we find useful to work with dimensionless momenta. The fact that we work with k (instead 
of with k) also implies that a factor kj^ ap pear s in front of every component of the Fourier components of the 
energy-momentum tensor. Concerning Eqs. (C.l) different comments are in order. First of all we can notice that 
at small scales (i.e. k > 1) the Fourier components of the energy momentum tensor are exponentially suppressed as 
exp [—k] = exp [—k/kg]. In other words the typical scale of the knot acts as ultraviolet cut-off. At large scales (i.e. 
for the infrared branch of the spectrum k < kg) it is of some interest to expand the obtained expressions and we find 
that 


Too(fc) 

Z'xx (^k) 
Tzz{k) 
Txz{k) 




TT + 1 


TT^Lik^y 2 4 


1 -kK-kO(K^) 


>ik) = 


Bl 


n^Likl V 2 5 


1 -I- K + 0{k^) 


Bl 


n'^L'tk^ 


— 3n^ 

• 

V 2 20 

1 -1- K -1- 0{k^) 


0(k), Tyz,{k) = . Oin), Tz;y{k) = 


Bl 


'n'^L^gkl 








(C.2) 


where O is the Landau symbol. We see that in the infrared limit (k < 1) the off-diagonal terms in the energy- 
momentum tensor are smaller than the diagonal ones. Of course at very small scales (smaller than the size of the 
knot ) the diagonal and off-diagonal terms are damped in the same way. We will use often the following notation 


T’oo(fc) = 

Txx{k) ~ 
Tzz{k) = 


Bl 


ir'^Ljk^ 


Tyy{k) = 


■K + 1 


Bl 


1 + k + 0{k^) 


fc3 V 2 5 


Bl 


■n ( h — 3n^ 


20 


1 -|- K -I- 0{k^) 


1 -I- K -I- 0{k^) 


(C.3) 


where we kept (in square brackets) the leading infra-red contribution and (outside the brackets) the leading (ultra¬ 
violet) damping term. The mathematical reason for this notation is very simple. In doing the ultra-violet expansion 
of the Fourier transformed components of Ty(fc) we can clearly see that the leading exponential damping factorises 
in the sense that each term appearing in the various components of is always multiplied by one of the different 
functions reported in Tab. I. The common feature of all these functions is the exponential damping which becomes 
effective at small scales and which is connected with the analytical form of the components of H which (in real space) 
all share the 1/(7?^ + 1)^ behaviour. Notice finally that, order by order in k the trace of vanishes. To lowest 
order, from Eqs. ( |C.3| ) one can immediately check that = 0 + 0{k). 




/3(k) = WT(l + «)e"" 

74 (k) = ^ ^ (3 + 3k + K.^)e ^ 

h{k) = ■^^J f (15 + 15 k -I- 6k^ -I- k')e 

h{k) = (105 -h 105 k -h 45 k^ -h 103 k'’ -h k^)e~'^ 


TABLE I. We report the coefficients appearing in the Fourier transform of the Cartesian components of the en¬ 
ergy-momentum tensor. We recall that, in our notations, k = |fc|/fcs. 
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